Bose-Einstein and Fermi-Dirac Interferometry 

in Particle Physics 



Gideon Alexander 



School of Physics and Astronomy 
Raymond and Beverly Sackler Faculty of Exact Sciences 
Tel-Aviv University, Tel-Aviv 69978, Israel 



Abstract 

The appUcation of the Bose-Einstein and Fermi-Dirac interferometry to multi-hadron 
final states of particle reactions is reviewed. The underlying theoretical concepts of parti- 
cle interferometry is presented where a special emphasis is given to the recently proposed 
Fermi-Dirac correlation analysis. The experimental tools used for the interferometry anal- 
yses and the interpretation of their results are discussed in some details. In particular 
the interpretation of the dimension r, as measured from the interferometry analysis, is 
investigated and compared to that measured in heavy-ion collisions. Finally the similarity 
between the dependence of r on the hadron mass and the interatomic separation on the 
atomic mass in Bose condensates is outlined. 
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1 Introduction 



The two-particle intensity interferometry {intensity correlation) method has been worked out 
and exploited for the first time in the 1950's by Hanbury-Brown and Twiss (HBT) [1-3] to cor- 
relate intensity of electro-magnetic radiation, arriving from extraterrestrial radio-wave sources, 
and thus measure the angular diameter of stars and other astronomical objects. This method 
rests on the fact that bosons obey the Bose-Einstein statistics so that the symmetrisation of 
the multi-particle wave function affects the measured many-particle coincident spectra to lead 
to an enhancement, relative to the single particle spectra, whenever they are emitted close by 
in time-phase space. Thus the interest in the Bose-Einstein correlation properties was not only 
restricted to their fundamental quantum theoretical aspects but above all as a tool in the study 
of the properties and spatial extent of sources emitting radiation and particles in the fields of 
astronomy, and non-perturbative QCD aspects of particles and nuclei interactions. 




Figure 1: The two indistinguishable diagrams that describe the emission of two identical bosons, tti 
and 7r2, emerging from the two points, ri and r2, which He within an emitter volume, and are detected 
at the positions xi and X2- 



The difference between the intensity interferometry and the conventional amplitude interfer- 
ometry can be illustrated [4,5] with the help of Fig. 1. Consider a finite source which emits two 
indistinguishable particles from the positions ri and r2 which are later observed at positions xi 
and X2- In an amplitude interferometry the positions xi and X2 could be a slit through which 
the emitted particles pass. The particles could then produce an interference pattern which 
will depend on the relative phase of the particle's amplitude as measured at xi and X2- In an 
intensity interferometry a normalised correlation function K (1,2) of the particles 1 and 2 is 
formed from the average number (ni^2) of counts where Ui and n2 are detected simultaneously 
at Xi and X2: 

The correlation function is thus proportional to the intensity of the particles at xi and X2- 
Because of the symmetrisation of their wave function, identical particles can have a nonzero 
correlation function even if the particles are otherwise noninteracting. 

As can be expected intensity interferometry is closely related to the amplitude interferometry 
which essentially measures the square of the amplitudes Ai and A2 falling on the detectors Xi 
and X2: 

+ = |Ai|2+|A2|' + KA2 + AiA^) , 

where the last term is the fringe visibility denoted by V, which is the part of the signal which 
is sensitive to the separation between the emission points. Averaged over random variations its 
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square is given by the product of the intensities landing on the two detectors: 

{V) = 2(|ylin^2|') + {AfAD + {AlAf) ^2{hh) . 

Since the last two terms of the (1/^) expression vary rapidly and average to zero, (K^) is pro- 
portional to the time averaged correlation of the product of the two intensities (see e.g. Ref. [6]). 

The Bose-Einstein correlation (BEC) analysis method was for the first time introduced in 
1959/60 by Goldhaber, Goldhaber, Lee and Pais (GGLP) [7,8] to the hadron sector in order 
to study the time-space structure of identical pions produced in particle interactions. Their 
observation of a BEC enhancement of pion-pairs emerging from pp annihilations at small rela- 
tive momenta was parametrised in terms of the finite spatial extension of the pp source and the 
finite localisation of the decay pions. The momentum range of the enhancement could then be 
related to the size of the particle source in space coordinates. Since then the boson interferome- 
try has been used extensively in variety of particle interactions and over a wide range of energies. 

In the mid 1990's the boson intensity interfcrometry measurement method has been ex- 
tended to the spin 1/2 baryon sector utilising the Fermi-Dirac statistical properties and in 
particular the Pauli exclusion principle [9] . Already the results from the first attempt to study 
the Fermi-Dirac correlation (FDC) of the AA and AA pairs, emerging from the decay of the 
gauge Z° boson, indicated that the spatial extension of the hadron source is strongly dependent 
on its mass. 




Figure 2: An early compilation of values obtained for the emitter dimension r^ms) extracted from 
interferometry analyses of identical charged pion-pairs produced in heavy ion collisions are displayed 
as a function of A^/^ where A is the atomic number of the projectile. The values are taken from 
Rcf. [10] and whenever more than one value is quoted for a given projectile nucleus, the plotted data 
points represent the average weighted values. The solid line represents the relation r rms — 

1.2 X A^/^ 

fm. 

In parallel to the BEC studies of hadronic final states produced in particle reactions, ex- 
tensive pion-pair correlation analyses were and are applied to heavy ion collisions in order to 
explore the characteristics of these rather complex reactions [6,11,12]. A compilation, taken 
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from Ref. [10], of the spatial dimension, rrms, extracted from the BEC of two identical charged 
pions emerging from nucleus-nucleus collisions, is shown in Fig. 2 as a function of A^^"^ where 
A is the atomic number of the projectile nucleus. As seen in Fig. 2, Vrms increases more or 
less linearly with A^^^. Moreover, as seen from the figure, the line Vrms = 1.2X74^/^ fm, which 
corresponds to the dependence of the nucleus radius on its atomic number A, follows rather 
well the BEC derived spatial dimensions. This observation, as well as the results from the HBT 
analyses in astronomy, led to the expectation that BEC should also render information on the 
geometrical extension of the hadron emitter source in particle collisions. 

The main purpose of the present report is to cover the field of BEC analyses applied to 
hadrons produced in particle collisions and only refer to heavy ion reactions whenever a com- 
parison between the two classes of reactions is of a particular interest. A special emphasis is 
given to recent developments in the particle interferometry such as the inclusion of identical 
fermion correlations and the new results concerning the spatial dimension dependence on the 
hadron mass. Furthermore, the report is aimed mainly to cover the experimental aspects of 
the BEC and FDC interferometry and their direct physics implication and achievements while 
trying to restrict the theoretical discussions to essentials so as not to overburden the non-expert 
reader. 

The organisation of this report is as follows. The basic concepts of bosons and fermions in- 
terferometry are discussed in sections 2 and 3. In section 4 details concerning the experimental 
aspects of the correlation measurements are described. Results obtained from the 1-dimensional 
correlation analyses are summarised in section 5. The observation that the spatial dimension 
deduced from the correlation measurements depends on the hadron mass is evaluated in section 
6 and conceivable grounds for this behaviour are discussed. Results obtained from the multi- 
dimensional correlation analyses are dealt with in section 7 and a relation between interatomic 
separation in Bose condensates and the spatial dimension of hadrons emerging in high energy 
reaction is illustrated. In section 8 the question concerning the interpretation of the BEC and 
FDC deduced geometrical dimension is discussed. Finally a summary is presented in section 9. 



2 Basic concepts in particle interferometry 
2.1 The Bose-Einstein correlation of two hadrons 

In deriving the main expressions used for the boson interferometry we follow closely Ref. [13]. 
As it is well known in quantum mechanics the interchange of two out of N indistinguishable 
bosons does not change the wave function describing this multi-boson state. This feature of 
the Bose-Einstein statistics means that the state ^ has the property that 

v[/(l,2,...,iV) = M/(2,l,..-,iV) , 

which leads to an interference term in that enhances near in time-phase space the produc- 
tion of indistinguishable bosons. Let us first consider a source of discrete emission points, pi, 
each characterised by a probability amplitude Fj (r) in the 3- vector phase space 

F,(r) = ^^(r-r,) . 
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Next we introduce the central assumption pertaining to the BEC effect namely, the chaotic 
or the total incoherence limit, which corresponds to the situation where the phases of the 
production amplitudes wildly fluctuate in every point of space. In this limit all the phases can 
be set to zero. If V'k(i') is the wave function of the emitted particle, then the total probability 
P(k) to observe the emission of one particle with a 3- momentum vector k is given by summing 
up the contributions from all the i points, that is 



P(k) = J2\p^ij{r, 
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For simplicity we will further use plane wave functions ■^k e**^'*'"''"''''^ where in the incoherent 
case we can set = 0. Next we replace the sum by an integral so that 

P(k) = I \p{r)\'£r. (2) 

The probability to observe two particles with momenta ki and k2 is 

P(ki,k2) = /iV'MlXrOlXr^jlVrA, (3) 

where ■0i,2 — "01,2 (ki, k2, ri, is the two-particle wave function. 

Taking incoherent plane waves, then for two identical bosons the symmetrised -01,2 is of the 
form 

J^s _ rgi(kiri+k2r2) , gi(kir2+k2ri)l 

•^1,2 ^ [ J \ J 

so that 

\ri,2\^ = l + cos[(ki-k2)(ri-r2)] = 1 + cos[Ak(ri - r2)] . (5) 
Using Eqs. (2), (3) and (5) one can define a second order correlation function 

r(l. ^^^'^^'^ - 1 ^ /cos[Ak(ri - r2)]|p(rO|Xr2)PdVidV2 ... 

= P(kOP(k2) - ' + P(kOP(k2) ' 

where Ak = ki — k2. Assuming that the emitter extension p(r) is localised in space and time 
then it follows that when Ak = the last term of Eq. (6) can vary between the values to 
1. Prom Eq. (6) one obtains after integration (Fourier transformation) 

C72(Ak)= 1 + |p(Ak)|2 . 

In many of the two-boson one dimensional BEC analyses one uses the Lorentz invariant param- 
eter Q, defined as = Ql = ~{Qi~Q2y = M2 — Here gi, 52 and M| are respectively 
the 4-momentum vectors and the invariant mass squared of the two identical bosons of mass /i. 
Thus one obtains 

C2{Q)^ 1 + |p(g)r- (7) 

Assuming further that the source is described by a spherical symmetric Gaussian density dis- 
tribution of emitting centres 

p(r) = p(0)e"^ , 
then the Bose-Einstein correlation function assumes the form 

C2(Ak) = 1 + e-'-o^'^' . 



In terms of the variable Q and the dimension rg, introduced by Goldhaber et al. [7,8], the 
correlation function in the completely chaotic limit is equal to 

C2{Q) = 1 + e-'-oQ' . (8) 

In the completely coherent case it can be shown [14] that C2{Q) = 1. In order to accommodate 
those cases where the source is not completely chaotic one introduces a chaoticity parameter 
A2 which can vary between the value 0, corresponding to a complete coherent case, to the value 
1 at the total chaotic limit. Thus Eq. (8) is transformed to the GGLP form 

C^iQ) = N{1 + X,e-'^oQ') , (9) 

where is added as a normalisation factor. Since the strength of the BEG effect depends 
also on the experimental data quality, like the purity of the identical boson sample, A2 is often 
also referred to as the BEG strength parameter. In the following, unless otherwise stated, 
we will denote by r the dimension values obtained from BEG analyses which used the GGLP 
parametrisation, that is r = tq- Two examples of a typical behaviour of the correlation function 
C2{Q) of identical charged pion-pairs are shown in Fig. 3. The first are the results of OPAL [15] 
where the pion-pairs were taken from the hadronic decays and the second reported by the 
ZEUS collaboration [16] in their study of the deep inelastic ep scattering produced at the HERA 
collider. 




Figure 3: The t:^t:^ BEC as a function of Q. Left: OPAL results [15] obtained in the hadronic 
decays. The sohd and dashed Hnes represent the fit results of Ref. [17] respectively with and without 
the inclusion of final state interactions (see Sec. 4.2). Right: ZEUS results obtained in deep inelastic 
ep scattering with momentum transfer of > 110 GeV^ [16]. The regions of Kg and which were 
omitted from the fit of the correlation function C2{Q) are indicated. 



2.2 The Kopylov-Podgoretskii parametrisation 

In addition to the Goldhaber parametrisation, given by Eq. (9), another favourite parametri- 
sation is that proposed by Kopylov and Podgoretskii (KP) [18-20] which corresponds to a 
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radiating sphere surface of radius tkp with incoherent point-hke oscillators of lifetime r, namely 



C{qt, go) = 1 + A[4j2(g,r^p)/(g,r|p]/[l + {qorf] , (10) 

where Ji is the first-order Bessel function. Here go is equal to \Ei — i?2| and g* = |q x k|/|k| 
where k = ki+k2, the sum of the two hadron momenta and q = ki— k2 is the difference between 
them. This correlation expression is not Lorentz invariant and its variables are calculated in 
the centre of mass system (CMS) of the final state hadrons. The correlation function C{qt,qo) 
can be written in terms of the Goldhaber geometrical parameter ro as [21] 

C{qt, go) = 1 + Aexp[-r^g2 - r^'Jij' - 1)] , (11) 

where 7 is the 7— factor of the identical boson pair. Prom a comparison between Eqs. (11) and 
(9) it is clear that the parameters and rxp have different interpretations. At small qt and 
go, however, the parametrisation (10) can be approximated to be 

Ciqt, go) = 1 + Aexp[-(r^p/2)2g2 - (q,rf] . (12) 

From Eqs. (11) and (12) one finds the approximate relation r^p — which is verified ex- 
perimentally. For example, in a study [22] of the BEC of the n'^'ri'^ system emerging from the 
decay of the Z° boson one obtained — 0.955 ± 0.012(6tat.) ± 0.015{syst.) fm as compared 
to rKP = 1.778 ± 0m3{stat.) ± 0.036(syst) fm. 



In contrast to the rather simple and straightforward BEC analysis offered by the GGLP 
parametrisation, the KP method expects a correlation enhancement at low, near zero, go values 
so that the analysis is carried out in several go energy slices. In addition the KP formalism 
has an extra free parameter, (gor)^, which has to be determined through a fit to the data. 
For these reasons the KP correlation investigations require higher statistical data samples than 
those needed for the BEC analysis in the GGLP method. 



2.3 Higher order Bose-Einstein correlations 



A Bose-Einstein correlation enhancement is also expected to be present in identical boson 
systems of more than two particles when they emerge from the interaction within a small time- 
space region. In the search for these so called, higher order BEC enhancements, one has to 
differentiate between those produced from the lower BEC order(s) and those who are genuine 
correlations. The normalised over-all inclusive correlations of n identical bosons is given by [23] 



Pn{Pl,P2 ■ ■■ Pn) 



^" Pl(Pl)Pl(P2) 



a 



da da 



Pl{Pn) 



— a 



dpidp2 ■ ■ ■ dpr, 



da 



^ dpi dp2 dpn . 



(13) 



where a is the total boson production cross section, pi{pi) and da/dpi are the single-boson den- 
sity in momentum space and the inclusive cross section, respectively. Similarly Pn{pi,P2 ■ ■ ■ Pn) 
and d'^ a / {dpidp2 ■ ■ ■ dpn) are respectively the density of the n-boson system and its inclusive 
cross section. The product of the independent one-particle densities Pi(pi)pi(p2) • • • Pi{Pn) is 
referred to as the reference density distribution, or reference sample, to which the measured 
correlations are compared to. Specifically the inclusive two-boson density ^2(^1,^2) can be 
written as: 

P2{PUP2) = Pi (Pl)Pl(P2) +^2(^1,^2) , (14) 
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where Pi(pi)pi(p2) represents the two independent boson momentum spectra and ^'2(^1,^2) 
describes the two-body correlations. In this simple case of two identical bosons the normalised 
density function i?2) defined by Eq. (13), already measures the genuine two-body correlations 
which here (see Sec. 2.1) is referred to as the C2 correlation function. Thus one has 

C2 = i?2 = (pi,P2) , (15) 

where 7^2 (P11P2) = -^2(pi,P2)/[pi(pi)pi(p2)] is the normalised two-body correlation term which 
in the GGLP parametrisation defined in Eq. (9), is equal to A2 exp(— Qirf). 

The inclusive correlation of three identical bosons, ^3(^1,^2, Pa), includes the three indepen- 
dent boson momentum spectra, the two-particle correlation K2 and the genuine three-particle 
correlation K^, namely: 

P3(Pl,P2,P3) = Pl(Pl)Pl(P2)pi(P3) +l]Pl(Pi)-f^2(Pj,Pfe) +-f^3(Pl,P2,P3) , (16) 

(3) 

where the summation is taken over all the three possible permutations. The normalised inclusive 
three-body density, is then given by 

Pl{Pl)Pl{P2)pl{P3) 

Here 



Pl(Pl)Pl(P2)pi(P3) ' ' Pl(Pl)Pl(P2)pi(P3) 

represent the mixed three-boson system in which only two of them are correlated and the three- 
boson correlation. In analogy to C2, one defines a correlation function C3 which measures the 
genuine three-boson correlation, by subtracting from the term which contains the two-boson 
correlation contribution. Thus 

C3 = R3-Rl,2 = 1+^3 (P1,P2,P3) , (18) 

which depends only on the genuine three-boson correlation. For the study of the three-boson 
correlation one often uses the variable Q3 which, analogous to the variable Ql, is defined as 

(3) 

where the summation is taken over all the three different i,j boson-pairs. Here M| is the 
invariant mass squared of the three-boson system and is the mass of the single boson. From 
the definition of this three-boson variable it is clear that as Q3 approaches zero so do all the three 
related qij values which eventually reach the region where the two-boson BEG enhancements 
are observed. It has been shown [24] that the genuine three-pion correlation function C3((53) 
can be parametrised by the expression 

C3(g3) = 1 + 2Xse-'isr's , (19) 
where A3 is the chaoticity parameter which may assume a value between zero and one. 
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The method outhned here for the extraction of the three-boson BEC enhancement can in 
principle be extended to higher orders, however it becomes too cumbersome to be of a practical 
use and on top of it, it requires very high statistics data. For these reasons two other approaches 
have been advocated and utilised experimentally. In the first, one measures experimentally the 
over-all correlation, as defined by i?„ in Eq. (13), and then with the help of various models 
one tries to extract the higher order genuine BEC parameters r„ and A„. Such an approach is 
adopted for example in references [25,26]. In the second approach one analyses the data in terms 
of correlation functions [27, 28] known as factorial cumulant moments or, in their integrated 
form are referred to as the semi-invariant cumulants of Thiele which constitutes an important 
element in proving the central limit theorem in statistics (see e.g. [29]). Having n particles 
in a given domain then the corresponding factorial cumulant moments will be different from 
zero only if a genuine n-particle correlation exists. The shortcoming of the cumulant approach 
is however the fact that its results for a given higher order {n > 3) pertain to the overall 
genuine correlation present in the data so that the partial contribution of the genuine BEC, 
parametrised by r„ and A„ is not transparent. 

2.4 Bose-Einstein correlation in two and three dimensions 

In the GGLP interferometry analysis the emitter shape is considered to be of a spherical shape 
with a Gaussian density distribution. A method to explore the possibility that the time-space 
extent of the particle emission region deviates from a sphere, and in fact is characterised by 
more than one dimension, has been recently proposed [11,30,31]. To this end the BEC analysis 
is carried out in the Longitudinal Centre-of-Mass System (LCMS) shown schematically in Fig. 
4. This coordinate system is defined for each pair of identical bosons as the system in which 




Figure 4: The Longitudinal Centre-of-Mass-System coordinates (taken from Ref. [32]). 

the sum of the boson-pair 3-momenta pi -|- p2, referred to as the out - axis, is perpendicular to 
the thrust (or jet) direction of the multi-hadron event defined as the z - axis (= long - axis). 
The momentum difference of the pion-pair Q is then resolved into the longitudinal direction 
Qz = Q\\ parallel to the thrust axis, to the Qout - axis which is collinear with the pair momenta 
sum and the third axis, Qside which is perpendicular to both and Qout- In this system the 
projections of the total momentum of the particle-pair onto the longitudinal and side directions 
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are equal to zero. The total value is given by 

Q' = Ql + Qlde + QUl-P'). where ^3 = ^if^^f^ . (20) 
Here the indices 1 and 2 refer to the first and second boson. Since pi^z — ~P2,z one has 

Qz = Q\\ = Pl,z - P2,z = '^Pi,z = ■ 

The difference in the emission time of the pions, which couples to the energy difference between 
them, appears only in the Qout direction. The 3-dimension correlation function is given by 

C2{Qz,Qs^cle,Qout) = 1 + A2e-(^^«^ + ''^^'^eQLe + r^^tQLt) . (21) 

In many cases due to lack of sufficient statistics, one wishes to reduce the number of pa- 
rameters to be fitted in the correlation function by defining the transverse component rr in the 
LCMS to be — r^^^g + r^„( corresponding to 

Qt — Qlut + Qlide ■ 

Thus the correlation function, which is fitted to the data, is of the form 

C2{Qz,Qt) = 1 + A2e-('-'«' + '•tQt) , (22) 

where r^, estimated from Eq. (22) as approaches zero, is the longitudinal geometrical radius 
and rT is composed of the transverse radius and the emission time difference. The experimental 

findings in heavy ion collisions [33,34], in hadron-hadron reactions [35,36] and e^e' annihila- 
tions [37-39], verify the theoretical expectations of the Lund string model [40,41] that the ratio 
tt/tz is significantly smaller than one (see Sec. 4.1.2). 

The assignment of a well defined physical direction in the LCMS, such as thrust axis, allows 
to study in the framework of the BEC analysis an additional meaningful variable namely, the 
average transverse mass ttit of the two identical hadrons. This transverse mass is defined as 

1 . , 1 r 



rriT = ^[mi,T + m2,T ] = 2iv"^i+^i'^ + V"^2+P2,tJ , (23) 
where the indices 1 and 2 refer to the first and second hadron. 



3 Fermi-Dirac correlation 

The application of two-boson correlation to the estimation of an r dimension has recently been 
extended to identical pairs of fcrmions (baryons) [9]. This extension is based on the Fermi-Dirac 
statistics feature which prohibit the total spin to have the value S = Si + S2 = 1 when the two 
identical fermions are in an s-wave {£ = 0) state. Thus in the so called Fermi-Dirac correlation 
(FDC) method one can, similarly to the BEC analysis, study the contributions of the S — and 
S — 1 states to the di-fermion system as Q approaches zero where, in the absence of di-baryon 
resonances, only the s-wave state survives. 

The estimation of r from the rate of depletion of the S = 1 population as Q ^ can be 
achieved in two ways. The first consists of a direct measurement of the relative contributions of 
S — 1 and S — spin states to the di-baryon system as a function of Q (see next section) . An 
alternative method consists of a measurement of the di-baryon density decrease as Q approaches 
zero, assuming its origin to be due to the Pauli exclusion principle. This second method is 
equivalent to that used in the BEC analysis of two identical bosons. 
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3.1 The spin-spin correlation 



If Fo{Q) and Fi{Q) are respectively the fraction of the S — and the S — 1 contributions, at 
a given Q value, to the di-baryon system then one can e.g. study two correlation functions, 
Cs=o{Q) and Cs=i{Q), defined by the ratios: 



CoiQ) - 



2Fo(g) 



and Ci(g) = 



2Fi(g)/3 



Fo{Q)+F^{Q)/3 Fo(g)+Fi(g)/3 ' 

where Fi{Q) is divided by 3 to offset the statistical 25 + 1 spin factor. At high Q values, where 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
Q (Arbitrary Units) 



§1.8 
ll.6 
1.4 
1.2 
1 
0.8 
0.6 
0.4 
0.2 




(b) 




1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 

y 



Figure 5: (a) Schematic view of the FDC functions Co and Ci dependence on Q, for the S = and 
S = 1 states, of an identical spin 1/2 two-baryon system. The dashed hnc represents the expectation 
for a baryon anti-baryon system in the absence of resonances, (b) The {l/N)dN/dy behaviour of the 
S = and S = 1 states of a two identical baryons sample. 



the highest angular momentum imax is large, one may expect that Cq = Ci = 1 corresponding 
to a statistical spin mixture ensemble. On the other hand when Q = 0, one will observe 
Co ~ 2 and Ci f« 0. Thus the Co behaviour as a function of Q is similar to the BEG function 
C2(g), of two identical mesons, which rises as the imax 0. As for the di-baryon correlation 
function Ci(g) there is no parallel case in the identical charged di-boson system. One further 
assumes the emitter to be a sphere with a Gaussian distribution so that, analogous to the BEG 
parametrisation, the di-baryon correlation functions can be parametrised as [9]: 

Co(g) = 1 + Xe-'"'^' and Ci(g) = 1 - Ae""'^' . 

Here the A parameter, which can vary between and 1, measures the strength of the effect and 
is mainly sensitive to the purity of the data sample. A schematic view of the dependence of 
Co and Ci on Q is shown in Fig. 5a. The study of either Cq or Ci can render a value for r 
however, since the diminishing contribution of the S = 1 state is responsible for the variation 
of the correlation function as Q decreases, one usually analyses Ci. Here it should be noted 
that unlike the BEG effect which can be extended to many-boson states, the FDG is limited to 
two identical fermions since already the third one is in an ^ > state (Pauli exclusion principle). 
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A general method for the direct measurement of the total spin composition of any two spin 
1/2 baryons system which decay weakly has been outlined in reference [9]. Here wc illustrate 
this method (later referred as Method 1) in its application to the AA and the AA systems. The 
simplest measurable case is the one where the A decays into pi + 7i~ and the other A, or A, to 
P2 + TT"*" where pi stands for the decay proton or anti-proton both of which we will further refer 
to as protons. Let us further denote by y the cosine of the angle (in space) between pi and 
P2 and (y) as its average, defined in a system reached after two transformations. The first, to 
the CMS of the di-A pair and the second where each proton is transformed to the CMS of its 
parent A. To note is that at Q = the second transformation is superfluous. 

In the decay of a single A, the proton angular distribution with respect to the spin direction 
in the centre of mass, is given by [42] : 

dw/dcos9 (X 1 — aACOsO , (24) 

where the parity violating parameter a\ = —aj was measured experimentally [43] to be 
0.642 ± 0.013. Setting a; = cos 6* the average {x\) is given by: 

(xa) — - — — / d(j) (1 — aAx)xdx where N — (1 — aAx)dx . (25) 

In particular, the average over the whole angular range from Xmin = — 1 to Xmax = +1 yields 

(xa) = -0.214 ± 0.004 and (x^) = +0.214 ± 0.004 . 

General arguments do not allow dN/dy distribution of a two spin- 1/2 hyperon system at 
threshold to have a y dependence higher than its first power. This means that dN/dy is of the 
form 

dN/dy ^A[l + Bs y] . 

The factor Bs can then be determined independent of the total angular momentum J value at 
Q = or nearby from the value of (y) using the Wigner-Eckart theorem [44,45]. Thus for the 
AA and AA pairs one obtains 

^5=0 = -9 (xa)^ = -0.4122 and ^5=1 = +3 (xa)^ - +0.1374 . 

As a result one then has for the AA(AA) system at, or very near, its threshold 

dN/dy\g^^ oc 1 - 0.4122y and dN/dy]^^^ oc 1 + 0.1374?/. (26) 

To note is that the £ = 0, S = 1 state is forbidden by the Pauli principle. These two distinctly 
different distributions, are shown in Fig. 5b. In a similar way one obtains for the AA system 

dN/dy\s^^ oc 1 + 0.4122y and dN/dyl^^^ oc 1 - 0.1374y . (27) 

Even though the Wigner-Eckart theorem is applicable to the di-baryon system at Q = 0, it is 
shown in Ref. [9] that Eqs. (26) and (27) can also be applied to at Q > as long as the protons 
emerging from the A decays are non-relativistic. That this method can be extended even to 
higher Q values has been pointed out in Ref. [46]. If the parameter e is defined as the fraction 
of the S = 1 contribution to the two-baryon system then it can be measured as a function of 
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Q by fitting the expression 



dN/dy = (1 - e) dN/dy\s=o + e dN/dy\s=i , (28) 

to the data. In the case of a statistical spin mixture, e is equal to 0.75 which yields a constant 
dN/dy distribution. This spin analysis method has the advantage that it directly measures 
the S =1 depletion as expected from the Pauli principle. In addition it does not require any 
reference sample which, as is well known from the BEC analyses, is the major contributor to 
the systematic errors of the measurement (see Sec. 4.1). The main disadvantage of the direct 
spin measurement is its need for a rather large data sample. In addition, it is limited to spin 
1/2 baryons which decay weakly into hadrons and thus is not applicable, for example, to pairs 
of protons or neutrons. 



3.2 The phase space density approach 

The probability to observe two particles, emitted from an interaction with 3-momenta ki and 
k2, is proportional to the square of its total wave function \I'i 2 which is equal to the orbital 
wave function ■0i 2 times the spin part. In the case of identical spin 1/2 baryons, the function 
^^1^2 should be anti-symmetric under the 1^2 exchange. Assuming plane waves and a 
completely incoherent emitter (i.e. the arbitrary phases can be set to zero) the symmetric and 
the anti-symmetric orbital wave functions are given by 



%/2 



gi(kiri+k2r2) _j_ gi(kir2+k2ri) 



and 



1 

71 



gi(kiri+k2r2) _ gi(kir2+k2ri) 



From these, as is shown in Sec. 2.1, one obtains for the symmetric orbital wave function, 
applicable to two identical hadrons 



whereas one obtains 



a |2 
2I 



1 + cos[(ki - k2)(ri - r2)] = 1 + cos(AkAr) , 



1 - cos[(ki - k2)(ri - Y2)] = 1 - cos(AkAr) , 



for the anti-symmetric orbital wave function. If we further consider a source with a spherical 
symmetric Gaussian density distribution [7,8] 



/(r) oc e-'-'/(^'-o) , 



then the rate will be 



s |2 
2I 



1 + e-''o 



and 



a |2 
2I 



1 - e '^o 



Because of the Fermi-Dirac statistics the symmetric orbital part of the di-baryon system, |V'i,2P) 
is coupled to the anti-symmetric spin part, that is = 0, whereas the anti-symmetric orbital 
part |'?/'i2p is coupled to the symmetric S = 1 spin part. If wc further assume that at high Q 
values we face a statistical spin mixture ensemble, where the probability to find a di-baryon 
in an S state is proportional to 25" -|- 1, we finally obtain for the emission rate, when properly 
normalised that 



1^1,: 



0.25[1 + e-'-o^'^' + 3(1 - e-'o^""^)] = 1 - 0.5e-'o 



r2Ak2 



2Ak2 
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or in its Lorentz invariant form 



^i,2(Q)|' = l-0.5e-'-'«' . 



Thus one expects at Q = a reduction in the correlation function to half of its value at the 
high Q region provided the following conditions are satisfied: 

• At high Q values one has a statistical spin mixture ensemble; 

• The di-baryon emitter is completely incoherent; 

• Absence of resonance states at low Q values; 

• Final states interactions can be neglected. 

To accommodate in this FDC analysis method (further referred to as Method II) the cases 
where the emitter is not completely incoherent one introduces, as in the BEC case, a chaoticity 
factor A, that can vary between and +1, so that the correlation function assumes the form 



where N represents a normalisation factor. 

4 Experimental procedure and data analysis 

Well suited reactions for a BEC analysis are those which lead to multi-hadron final states where 
the correlations due to resonances and conservation laws, like energy-momentum and charge 
balance, have minor effects. Since in these reactions the fraction of pions is the highest one, in 
many of the analyses one assumes that all the outgoing particles are pions and the contamina- 
tion from other hadrons are accounted for by a proper correction factor and/or by an increase 
of the systematic error. The correlation study of kaon and proton pairs require special hadron 
selection criteria which tend to reduce the data statistics and introduce larger systematic errors. 

In the BEC analysis the space density of the data hadron pairs dependence on Q is compared 
to a reference sample distribution which serves as a yardstick. The correlation function, of the 
form given in Eq. (9), is then fitted to the ratio of the data density distribution to that of 
the reference sample. To account in the GGLP parametrisation, for the so called 'long range 
correlations' due e.g. to energy and momentum conservation, Eq. (9) is often modified to 
include linear and quadratic terms in Q, namely 



where 5 and r] are free parameters to be determined by the fit to the data. In assessing the 
fitting results obtained for the r and A parameters one has to be aware of the fact that in 
many cases these are not independent and a correlation between the two does exist. The weak 
part of the BEC analysis is undoubtedly the fact that it depends rather strongly on the chosen 
reference sample which thus is the main contributer to the over all systematic errors associated 
with the fitted parameters. The different available reference samples and the effect of the final 
state interactions on the BEC analysis results are discussed in some details in the following 
sections. 



C{Q) = N{1 - O.SAe"'-'^') , 



(29) 



C{Q) = N{1 + Ae-'-''?')(l + 5Q + riQ^) , 



(30) 
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4.1 Choice of the reference sample 

As mentioned above the study of the second and higher order particles' BEC enhancements 
requires a yardstick against which they can be detected and measured as is also evident from the 
C2(kik2) definition given by Eq. (6). This yardstick is given by the so called reference sample 
that should be identical to the analysed data in all its aspects but free from Bose-Einstein 
or Fermi-Dirac statistics effects. An ideal solution to this requirement really does not exist 
but one has a choice of possibilities which satisfy approximately this requirement. These are 
divided mainly into two categories. Reference samples constructed out of the data themselves 
and those supplied by Monte Carlo generated samples which are subject to a full simulation of 
the experimental setup with its particle detection capabilities. 

4.1.1 Reference samples derived from the data 

The methods where the reference samples are constructed from the data themselves are often 
preferred as they are expected to retain many of the kinematic and dynamical data correla- 
tions such as those originating from charge, momentum and energy conservation as well as 
correlations arising from hadronic resonances. In the following we will briefly describe some of 
the data derived reference samples which, were and still are, utilised in the BEC analyses of 
pion-pairs. 

a) In many studies of the BEC of identical charged pion-pairs (tt^tt^) the simplest refer- 
ence sample is used. Namely, the one which is constructed out of the correlation of opposite 
charged pion-pairs (tt'^tt^) present in the same data sample. This choice however has the fol- 
lowing rather severe drawback. Whereas the tt^tt^ is a so called exotic system void of bosonic 
resonances, the origin of the tt^tt^ pairs may also be resonances where among them the most 
dominant one is the p(770). To overcome this deflciency the expression of the correlation func- 
tion, like that given by Eq. (9), is fitted to the experimental results only in the Q regions where 
the data is known to be free of resonances. 

b) Another frequently used reference sample is the one known under the name 'mixed-event' 
sample. In this method one couples two identical pions each originating from a different data 
event. In this way one is guaranteed that no BEC effects will exist in the sample but at the price 
that all other kinds of correlations, like those arising from kinematic conservation laws, are also 
eliminated. Furthermore, as long as the data analysed is produced at low energies where the 
the particles emerge to a good approximation isotropically, this event mixing procedure may be 
satisfactory. However at higher energies where the particles emerge in hadron-jets, the mixed 
event technique can only be applied to two-jets events where the mixing takes place between 
two events with a thrust (sphericity) axes lying very nearby or alternatively after one event is 
rotated so that the two events axes coincide. 

c) A third method in use for the generation of a reference sample, which avoids the need to 
rotate the event, is constructed by folding each data event along its sphericity or thrust axis so 
that the emerging hadrons are divided into two hemispheres. The entries to the reference sam- 
ple are then all possible identical pion-pairs belonging to different hemispheres. This method, 
as the former one, can only be applied to two-hadron jet events. 

Finally to note is that the choice of data derived reference samples for a BEC analysis 
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of kaon-pairs or for a FDC Method II analysis of baryon-pairs is much more restricted. For 

example, in the BEC analyses of K^K^ or K^Kg pairs the data K^K"^ pairs cannot be used 
due to the strong presence of the 0(1020) — * K^K^ decay which lies at the very low Q range 
where the BEC interference is near its maximum. 

4.1.2 Monte Carlo generated reference samples 

Modelling of the hadron production in particle reaction plays a central role in any experimental 
data analysis in high energy physics. Its aim is mainly the evaluation of the various experi- 
mental deficiencies due to the imperfection of the detection system such as the geometrical 
acceptance, angular resolution and separation of tracks as well as the limited particle identifi- 
cation capabilities. In addition, the hadron production simulation serves as a yardstick against 
which various physics phenomena and hypotheses can be detected and measured. In this last 
capacity modelling of hadron production is also used extensively in the study of particle corre- 
lations. 




electro-weak 



Figure 6: A schematic diagram for the production of hadrons in e^e annihilation. 



To illustrate the hadronisation process we show in Fig. 6 a schematic diagram for the pro- 
duction of hadrons in the process e^e^ /'-y* qq. The gg, which is not directly detected, 
results in a relative large number of hadrons. This final step process is governed by the strong 
interactions and can be divided into two parts. Firstly a parton cascade develops in which the 
q and q pair radiates gluons, which in turn may radiate additional gluons and split into new 
qq pairs. This process is followed by the hadronisation stage which emits the final observed 
particles. These two last stages arc described in the standard model of particles and fields in 
terms of the Quantum Chromo Dynamics (QCD) theory where the force field between partons 
is the so called colour field. 



At short distances and over short times the quarks and gluons can be considered as free 
particles where the perturbative QCD can be applied. This is not the case at the hadronisation 
stage where there is a need to rely on non-perturbative QCD models. One of the leading suc- 
cessful models, which currently is applied to many of the high energy particle reactions leading 
to hadronic final states, is the Lund string fragmentation model [47], shown schematically in 
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Fig. 7. In this Lund model the probabihty \M{qq — > to emit n-hadrons is propor- 

tional e~^^ where A is the colour field area spanned in time-space by the primary qq pair and 
6 is a real positive constant. 

In the simplest case where the q and q emerge in opposite directions, the colour field spanned 
between them is approximated by a massless relativistic string with constant energy density. 
These quarks oscillate back and forth and if given by collision enough energy the string will 
tear and a new pair of quarks will be produced. This process repeats itself until only ordinary 
hadrons will remain. In this model the produced hadrons are given a small transverse momen- 
tum with respect to the string axis. To note is that the production of quark pairs with nig > 
costs energy which is taken from the constant energy density string which means that they 
cannot be produced at the same location but will be separated by some distance proportional 
to their mass. The success of the Lund model stems not only from its ability to describe many 
of the high energy interaction features but also from the fact that it can be formulated stochas- 
tically as an iterative process and therefore is well suited for computer Monte Carlo simulation 
programs like the JETSET and its more recent advanced versions. 




Figure 7: The hadronic decay of a Lund model string spanning the time-space area A. 

Another approach to the hadronisation has been taken by the cluster model [48]. This 
model is based on the observation that partons generated in a branching process tend to be 
arranged in colour singlet clusters formed from qq pairs with limited extension in coordinate 
and momentum space. These clusters have masses typically of the energy scale at which the 
parton shower terminates. Very massive clusters decay first into lighter clusters pairs and then 
each cluster decays isotropically into observable hadron pairs with branching ratio determined 
by the density of states. This model is implemented in the HERWIG Monte Carlo program 
ready to be used in the experimental analyses. 

At this point it is important to note that although outgoing particle tracks generated in the 
Monte Carlo programs do have a charge ascription their Coulomb effects are not included. In 
addition, all generated particles are spinless. In the Lund model derived Monte Carlo versions 
a facility has been added so that the Bose-Einstein statistics effects present in identical boson 
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system can be incorporated [49]. Another deficiency of the Monte Carlo programs is the need 
to adjust quite a few free parameters, hkc the ratio between vector and scalar resonances and 
their production cross sections, to the experimental findings. As a consequence, one should 
proceed with caution when comparing a particular effect seen in the data with the Monte Carlo 
program prediction. Moreover, the setting of the many free parameters may well interfere with 
the possibility to test the validity of the assumptions and underlying building elements of the 
particular hadronisation model. 



4.2 Final state interactions 

In the BEC and FDC analyses, final state interactions (FSI) may also play a role. There are 
two major FSI types. The first is the Coulomb interaction which affects the charged hadron 
systems and the second, the strong interaction final state which is present in both charged 
and neutral hadron systems. So far many of the reported BEC analyses have included the 
Coulomb effect whereas the strong final state interactions have in general been avoided due to 
their complexity and the realisation that their effect on the A and r parameters is relatively 
small and can be absorbed in the overall systematic errors (see Sec. 4.2.2). 



4.2.1 The Coulomb effect 

In the Bose- Einstein correlation of identical charged hadrons the Coulomb repulsive force tends 
to reduce the enhancement signal. If C™^"*((5) denotes the measured correlation of two hadrons 
then its relation to the true correlation C2{Q) is given by 

Cr^^g) = C^iQ) X G2{Q) , (31) 
where G'2.{Q) is determined through the Gamow penetration factor [50] 

G2{Q) = ^,^J''l ^ where r^,,, = t^t^-^ . (32) 

Here are the charges in positron units of the hadrons, m their mass and a is the fine- 
structure constant. As seen from Eq. (32), the Coulomb effect increases as Q approaches zero. 
The measured correlation is corrected by the factor XjGiiff) which may produce an exaggerated 
BEC signal [51] and therefore proper caution has to be exercised in its application. A somewhat 
different Coulomb correction method has been advocated in [52]. To note is that the Coulomb 
correction for a di-pion system is rather small (see Fig. 8) and even at Q2 = 0.2 GeV does 
not amount to more than 2% and therefore in many reported BEC studies this correction was 
ignored. 



The Coulomb effect on the BEC analysis of three identical charged hadrons can be expressed 
[54] , to a good approximation, in terms of two-hadron Coulomb factors as 

^3(^3) = G2{Ql,2) X ^2(^1,3) X G2(Q2,3) ■ (33) 

Further improvements to Eq. (33) have been proposed (see e.g. Refs. [55,56]). Unhke the 
di-pion case the Coulomb correction for three equally charged pions at Qs — 0.25 GeV is not 
negligible and, as seen in Fig. 8, amounts to about 7%. 
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Figure 8: The Coulomb correction to the BEC correlation function [53]. (a) For two identical charged- 
pion systems as a function of Q2 and (b) for the three charged pion systems as a function of Q^. 



4.2.2 Strong final state interactions 

The simultaneous effect of both the Coulomb and the strong interaction scattering of two iden- 
tical charged pions on the BEC analyses have been lately worked out and reported in Ref. [17] 
where references to earlier studies are also included. The FSI of the strong type are limited, 
due to the short range of strong interaction, to the s-wave alone. For the vr^vr^ system the 
FSI dependence on Q is given by the well measured /^,r = 2 phase shift 5g {Q) which can be 
incorporated into the BEC function. 

An example for the inclusion of FSI in a BEC analysis Ref. [17] is illustrated in Fig. 3. In 
this figure the data points shown are the vr^vr''^ correlation function versus Q as measured by 
the OPAL collaboration [15] in the hadronic decays. The solid line represents the fit results 
to the BEC function including Coulomb and strong FSI whereas the dashed line is the outcome 
of a fit where the FSI were ignored. In a systematic study of the FSI in several BEC analyses 
applied to e^e~ pions data collected at energies on the Z° mass and below, it was found 
that in general the inclusion of the FSI tends to increase A and decrease the r value [17]. 



5 The 1-Dimensional correlation results 
5.1 The TTTT Bose-Einstein correlation 

The major information concerning the BEC features is coming from the vr^vr^ system which 
has been analysed in a large variety of particle reactions and over a wide range of centre of 
mass energies. In Table 1 results obtained from the 1-dimensional BEC analyses of two iden- 
tical pions produced in e^e~ annihilation in the energy range from 29 to 91 GeV are listed. 
The values for r2 and A2 are divided into two groups according to the reference sample type 
that was used. In Method I the reference sample chosen was the data vr^vr^ system whereas 
in Method II the reference sample used was either a generated Monte Carlo sample plus a full 
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Table 1: The two-pion GGLP emitter dimension r2 and the chaoticity parameter A2 deduced from 
BEC studies of e+e~ annihilations using one or two choices for the reference sample. In Method I 
the reference sample was the 7r"'"7r~ data sample. In Method II the reference sample used was either 
Monte Carlo generated events or a sample created by the event mixing technique. The given errors are 
the statistical and systematic errors added in quadrature. The values given are without a Coulomb 
correction. 



7r±7r± BEC 


Method I 


Method II 


Experiment 


^ [GeV] 








Ao 

/A2 


MARK II [24] 


29 


0.75 ± 0.05 


0.28 ±0.04 


0.97 ± 0.11 


0.27 ± 0.04 


TPC [57] 


29 






0.65 ± 0.06 


0.50 ± 0.04 


TASSO [58] 


34 


0.82 ± 0.07 


0.35 ± 0.03 






AMY [59] 


58 


0.73 ± 0.21 


0.47 ± 0.07 


0.58 ± 0.06 


0.39 ± 0.05 


ALEPH [60] 


91 


0.82 ± 0.04 


0.48 ± 0.03 


0.52 ± 0.02 


0.30 ± 0.01 


DELPHI [61] 


91 


0.83 ±0.03 


0.31 ± 0.02 


0.47 ± 0.03 


0.24 ± 0.02 


L3 [62] 


91 






0.46 ± 0.02 


0.29 ± 0.03 


OPAL [22] 


91 


0.96 ± 0.02 


0.67 ± 0.03 


0.79 ± 0.02 


0.58 ± 0.01 


7r°7r° BEC 


r2 [fm] 


A2 


r2 [fm] 


A2 


L3 [62,63] 


91 






0.31 ± 0.10 


0.16 ± 0.09 


OPAL [64] 


91 






0.59 ± 0.11 


0.55 ± 0.15 



detector simulation or a sample obtained by the event mixing technique. 

From the n'^Tr'^ BEC analyses results shown in Table 1 one can make the following observa- 
tions. Firstly the values of r2 are in general found to be higher in the Alcthod I analysis than 
in the Method II analysis whereas the A2 values are roughly the same. Secondly there is no 
clear evidence for an r dependence on the centre of mass energy. To note is also the fact that 
the values of r2 are smaller than 1 fm and their average in Method I is in the vicinity of 0.8 
fm. In Method II the r values arc less stable and fluctuate from experiment to experiment. It 
is generally believed that these fluctuations of the r values are just a reflection of the fact that 
different experiments used different BEC analysis procedures, namely in their data selection 
criteria and in their choice of the reference sample. This situation is also clearly seen in Fig. 
9 where r2 is plotted against the e'^e" centre of mass energy, Ecm- In particular to note are 
the three r2 values obtained in Method I for the tt^tt^ pairs emitted from the gauge boson 
where the OPAL value differs by several standard deviations from those reported by ALEPH 
and DELPHI. The L3 collaboration has not given an r value in the framework of the Method 
I analysis. Unlike the abundant results on the n'^n^ correlations the information on the BEC 
analyses of tt^tt^ pairs is very limited due to the experimental difficulties in identifying the 
two neutral bosons. At the same time there exists an interest in the two neutral pion sys- 
tem both in the framework of the string model and in the so called generalised BEC where 
isospin invariance is incorporated into the analysis. The two r values, obtained from the BEC 
analysis of the tt'^tt^ system, present in the hadronic Z° decay and listed in Table 1, do not 
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Figure 9: The emitter dimension r^±^± obtained from BEC analyses of e^e~ annihilation into hadrons 
versus the centre of mass energy of the e~^e~ system. The values were taken from references [22, 24, 
57-62] where all but TPC and L3 used a data reference sample (see also Table 1). 



allow presently a meaningful comparison to the string model and isospin invariance predictions. 

A compilation of n^n'^ BEC parameters obtained from analyses carried out in a variety of 
particle reactions, other than e~^e~ annihilation, is given in Table 2. Again no clear evidence 
can been seen for an r2 dependence on the type of reaction and/or energy. The r values re- 
ported are, like in e~^e~ annihilations, again near, or somewhat lower, than 1 fm. 

The BEC of the tt^tt^ pairs produced in the reaction e~^e~ W^W~ —>■ hadrons, where 
the pions are emerging from different W bosons, have been recently investigated in view of 
the concern that they might affect the W gauge boson mass estimation [75]. Whereas the 
preliminary results of DELPHI found some evidence for a BEC enhancement [76] , the L3 and 
ALEPH collaborations have not seen any effect of this kind [77, 78] . 

5.2 The K^K^ and /C^iTj systems 

The results from BEC analyses of the K^K^ system are meager in comparison to those re- 
ported from the tttt BEC analyses. This is mainly due to the fact that the rate of the charged 
kaons is much smaller than that of the charged pions which even in the hadronic Z° decays 
amounts only to ~ 2.4/17. An additional reason for the relatively poor kaon-pairs statistics 
stems from the need to remove from their sample contributions from protons, anti-protons and 
charged pions. This is achieved by applying hadron identification criteria which in most cases 
are effective only in a limited range of momentum. In the e~^e~ annihilations on the mass 
energy two K'^K^ BEC analyses have been reported [79,80]. Their extracted rK±K± values 
are listed in Table 4. 
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Table 2: A representative sample of the two-pion emitter dimension r and the chaoticity parameter A2 
deduced from BEC studies of hadron-hadron, lepton-hadron and 77 reactions leading to multi-hadron 
final states. The data selection and reference samples, as well as the application of the Coulomb 
correction, differ from experiment to experiment. In most cases the errors given are only the statistical 
ones. The systematic errors are typically of the order of 10 to 20% of the measured parameter values. 



n^n^ BEC analyses 


Parameter 


Reaction 


EcM [GoV] 


1-2 [fill] 


A2 


77 ^ /i [24] 
77 — > Btt^ [65] 


< 5 > 
1.6 - 7.5 


1.05 ± 0.08 
0.54 ± 0.22 


1.20 ± 0.13 
0.59 ± 0.20 


u{P)N h [66] 


8-64 


0.64 ± 0.16 


0.46 ± 0.16 


jip ^ h [67] 


23 


0.65 ± 0.03 


0.80 ± 0.07 


A +7; ^ // [(38] 


21.7 


0.83 ± 0.0(3 


0.33 ± 0.02 


pp ^ h [69] 
pp ^ h [70] 
pp ^ h [71] 


26 
27.4 
63 


1.02 ± 0.20 
1.20 ± 0.03 
0.82 ± 0.05 


0.32 ± 0.08 
0.44 ± 0.01 
0.40 ± 0.03 


pp^h [72] 

PP ^ h [73] 


1.88 

200 - 900 


1.04 ± 0.01 

0.73 ± 0.03 


1.96 ± 0.03 

0.25 ± 0.02 


e p ^ e h [74] 
e p ^ e h [16] 


6 < Q2 < 100 [GeV2] 
110 [GeV^] 


0.68 ± 0.06 
0.67 ± 0.04 


0.52 ± 0.20 
0.43 ± 0.09 



Unlike the K^K^ pairs, the origin of the K^K^ system is not only fed from strangeness 
±2 states but also from the boson-antiboson K^K (strangeness = 0) system which in turn 
may be fed from the decay of resonances. Neglecting CP violation^ we show in the following 
that a Bose- Einstein like threshold enhancement is nevertheless expected in the Q{Kg, Kg) 
distribution even when its origin is the K^T{^ system [81-84]. 



In general a boson-antiboson, BB, pair is an eigenstate of the charge conjugation operator 
C. In the absence of outside constraints, one can write the probability amplitude for a given 
charge conjugation eigenvalue Cn as follows: 



1 

72 



B:B 



Cn=±l 



^ B{p);B{-p)) ± 1 \B{py,B{-p 



(34) 



where p is the three momentum vector defined in the BB centre of mass system. In the Q — 
limit, where the BEC effect should be maximal, p is equal to zero and Eq. (34) reads 
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B:B 



c„=±i 



1 
2 



5(0); 5(0)) ± I 5(0); 5(0) 



(35) 



This means that, at Q = 0, the probability amplitude for the C„ = — 1 state (odd i values) is 
zero whereas the C„ = +1 state (even i values) is maximal. Here it is important to note that 

^The inclusion of CP violation effect in the following discussion is straightforward but with negligible conse- 
quences at the current experimental precision level. 
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these equations do hold for any spinless boson-antiboson pair such the K^K , the K'^K" and 
the TT+TT" systems. What does distinguish the K'^K from the other spinless boson-antiboson 
pairs is the simplicity by which one is able to project out the C„ = +1 and the C„ = — 1 parts 
of the probability amplitude. 

As is well known, the K" and the K mesons are described in terms of the two CP eigen- 
states, Kg with CPn = +1 and with CPn = -1. From this follows that, when the K^K° 
pair is detected through their Kg and decays, a definite eigenvalue C„ of the K'^K sys- 
tem can be selected. Thus, as Q approaches zero, an enhancement should be observed in the 
probability to detect Kg Kg pairs and/or K^K^ pairs (C„ = +1), whereas a decrease should 
be seen in the probability to find KgKl pairs (C„ = —1). The BEC can thus be analysed by 
forming for the eigenstate C„ = +1 and C„ = —1 respectively the correlation functions: 



P{K^Tf) P{K^Tf) 
where P stands for the production rate of a given two-boson state as a function of Q. 

These Cj^{Q) and C-{Q) dependence on Q, for A = 1. are drawn schematically in Fig. 10. 
As seen, when Q approaches zero, C{Q) splits into two branches. The first rises up to the vahie 
two and the other decreases to zero in such a way that their sum remains constant and equal 
to one at all Q values. This means, that if all the decay modes of the K^K pairs are detected 
and used simultaneously in the same correlation analysis then, according to Eqs. (34) and (35), 
no BEC effect will be observed at Q = 0. This however should not come as a surprise if one 
recalls that the K^K system is at the very end not composed of two identical bosons. 




Q (Arbitrary Units) 

Figure 10: Schematic behaviour of the correlation functions C{Q) of the system given by 

Eq. (36). The C„ = +1 probability ampUtude reaches the value 2 at the limit of Q = whereas 
the C„ = — 1 probability amplitude reaches the value 0. The sum of Cn = +1 and C„ = — 1 remains 
constant down to Q = 0. 



The interpretation of a BEC enhancement in the KgKg system is not straightforward as 

23 



its origin may also be the decay of the scalar /o(980) resonance. This resonance, which is 
mainly observed in its decay to pion-pair in the decay of the gauge boson [85], lies within its 
width at the K^K^ threshold of 995.3 MeV. Thus the low Q enhancement seen in the K^^ 
correlation at Q ~ may be formed from the Bose- Einstein effect and/or from the /o(980) 
resonance decay. A possible method to disentangle the two contributions is given by isospin 
considerations [86] which is outlined in the following section. Values for A and r obtained for 
the Kg system, produced in e+e~ annihilation on the Z° mass, are given in Table 4. 

5.3 Isospin invariance and generalised Bose-Einstein correlation 

In the sector of strong interacting hadrons where BEC and FDC analyses are performed, the 
isospin is a conserved quantity. In analogue to the so called generalised Pauli principle, which 
is the extension of that principle from two identical nucleons, the proton-proton and neutron- 
neutron pairs, to the proton-neutron system, one can generalise the Bose-Einstein statistics by 
the inclusion of the isospin invariance. This extension may relate for example, BEC results 
of the T^^T^^ pairs to those obtained for the 7r+7r° system. This extension was pointed out by 
several authors (see e.g. references [13,87-89]) and recently was worked out in details [86] to 
produce concrete relations between the BEC effects of hadron systems which are connected 
by isospin and are proposed to be tested experimentally. So far however the concept of the 
generalised BEC has not been verified as it involves pion systems which include one or more 
neutral pions which experimentally are hard to detect and identify. However if the gener- 
alised BEC assumption is confirmed then its effect on particles correlations is not negligible. In 
the following we briefly discuss some aspects of the generalised BEC features and consequences. 

First to note is that in many cases boson pairs, like the KK and tttt systems, are produced 
together with other hadrons, here denoted by X, from an initial state which is isoscalar to a 
very good approximation. These include e.g. pairs produced from an initial multi-gluon state, 
pairs from hadronic decays of J/ip and T, or pairs produced by decays. In some of these 
cases the initial state is however not pure 1=0 but has also some contamination of an 1=1 
component which is mixed in like in those processes where the and T decay to hadrons 
via one photon annihilation. According to the specific case methods can be applied to reduce, 
or even eliminate, this contamination. For example, the sTibsamplc of the C„ = — 1 quarkonia 
which decays into an odd number of pions assures, due to G-parity, that the hadronic final 
state is in an 1=0 state. Forming such a subsample from the hadronic Z^ decays is not useful, 
but on the other hand multi-hadron final states which originate from Z° decays into ss, cc and 
U) quark pairs, are in an 1=0 state. 

Several specific relations in the framework of the generalised Bose-Einstein statistics be- 
tween charge and neutral K-mesons have been proposed in reference [86]. Among others, of 
a particular interest are those relations which have a bearing on the analysis of the enigmatic 
scalar /o(980) resonance, whose mass and width he in the KK threshold region, and which has 
a long history regarding its nature and decay modes [43, 90-92] . For details concerning these 
relations the reader is advised to turn to reference [86]. 

Next we turn to the effect of the generalised BEC assumption on the tttt system. In this 
system there are more states and more isospin amplitudes than the two, 1=0 and 1=1 of the 
KK system. In the low Q region where only the s-wave contributes, the tt'^tt'*', tt'^tt" and 7r°7r° 
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amplitudes depend upon only two isospin amplitudes, 1=0 and 1—2. Their intensities in this 
region satisfy a triangular inequality [86] 



Y: V(2/3) • P% ^ (7r0^o)X] - ^(1/3) • P[i.„ ^ {n+7r-),X] 

X 

< E \/P[io - (7r±7r±)X] = Y: ^JP[io - (7r±7rO),X] < 

X X 

<Y V(2/3) • P[io - (7r%o)X] + ^(1/3) • P[z„ ^ (7r+7r-)eX] 



< 



(37) 



where the notation (7r7r)e is used to indicate that only even partial waves for the (tttt) final state 
are included in the sum. Here ig denotes an initial 1=0 state and X stands for the hadrons 
produced in association with the two earmarked pion-pair. 

In this last relation the middle line is of a particular interest as it states that at low Q 
values, where only the s-wavc survives, the BEC enhancement in the ti'^ti^ system should be 
equal to that present in the n^n^ pairs. Finally the BEC features of the tt^tt" system should 
not be far from that of the n'^n'^ state as long as the contribution of the tt^tt^ to the triangle 
inequality relation remains small. 



5.4 Observation of higher order Bose-Einstein correlations 

As mentioned in Sec. 2.3, the directly measured higher order BEC are essentially hmited to the 
case of three identical charged pions. In these analyses, the contribution from the two-hadron 
BEC was subtracted in each experiment with somewhat different method. The genuine tt^it^tt^ 
BEC were also measured in the reaction e~^e^ followed by an hadronic decay [93-95]. 

The experimental normalised distribution, which is marked by (a) in Fig. 11, is that found by 
OPAL [94] to describe the contribution to R-i{Qz)., as defined in Eq. (17), from the over-all 
three identical charged pion BEC enhancement. In the same figure the three-pion distribution, 
denoted by (b), corresponds to the BEC enhancement due to the lower order BEC correlations 
i.e. the two-pion system. A clear evidence is seen for the genuine three-pion BEC enhancement 
which can be extracted by subtracting the distribution (b) from that of (a) which then allows 
the evaluation of r-^. 

Results for the three-boson emitter dimension are listed in Table 3 where the values of rs 
and the ratio ro^jr2 arc given. As can be seen, it does not seem to exist a significant dependence 
of the r3 and r-ilr^ values on the type and/or energy of the particles reaction. 

In Ref. [24] a relation is derived between the two-pion emitter radius r-i and the three-pion 
emitter radius ra as extracted from a fit of RsiQs) to the over-all (non-genuine) three-boson 
correlation distribution. This relation, which is based on the Fourier transform of a Gaussian 
shape hadron source, is given by the inequality 

ra/VS < ra < rs/v^ . (38) 

so that r3/r2, lies in the range of 0.55 to 0.71. The inequality relation (38) is reduced to the 
equality 

rs = r2/V2 , (39) 
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Figure 1 1 : BEC results from an OPAL study of three identical charged pions detected in the hadronic 
Z° decays [94] . The over-all normalised BEC enhancement is described by the (a) distribution whereas 
distribution (b) corresponds to the three-pion BEC enhancement generated from the lower order BEC 
enhancement of the two-pion system. 



when the emitter radius is determined from the genuine BEC distribution parametrised by 
C'alQa) so that r3/r2 should exactly be equal to 0.71. The experimental results, shown in Table 
3, are seen to fulfil the expected relations between r2 and and in particular the results from 
the genuine three-pion correlation obtained from the hadronic Z° decay are consistent within 
one to two standard deviations with the expectation of Eq. (39). 



5.5 Experimental results from Fermi-Dirac correlation analyses 

Until recently FDC analyses of fermion-pairs were prohibited due to the low production rate 
of identical di-baryons in particle reactions at low and intermediate energies. This situation 
changed with the commission of the LEP collider, operating on the Z° mass energy, where 
each of its four experiments have accumulated some 3 to 4 million hadronic Z° decays. The 
inclusive decay branching ratio of ~0.34 A per event was sufficiently high to allow a FDC 
analysis of the A A(A A) and A A pairs which are free of Coulomb effects and are relatively easy 
identified. Results from these analyses have been reported by the OPAL [98], DELPHI [99] and 
ALEPH [100] collaborations using the spin-spin correlation method (see Fig. 12). The ALEPH 
collaboration has further repeated its analysis with Method II, i.e. with the phase space density 
approach and found the results of both methods to be consistent within errors (see Fig. 13). 
In the same figure is also shown the OPAL [101] recent p p FDC measurement carried out in 
the hadronic Z° decays. 

A summary of the measured di-baryon r values is given in Table 4 together with the average 
LEPl value for the pairs and values for the K"^ and Kg K'^ systems obtained from 

BEC analyses. To note is the fact, as already discussed in Sec. 5.2, that the pairs 
may also be the decay product of resonances like the /o(980), whereas the other hadron pairs 
listed in the Table form, so called, exotic states, that is they are in isospin 2 or strangeness ±2 
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Table 3: Results for and r3/r2 values obtained from BEC analyses of two and three identical 
charged pions. Values for genuine three-pion BEC are marked by (g) otherwise they are marked by 
(n). 



27r and Stt BEC analyses 



Reaction 



pp ^ h 
pp ^ h 
pp ^ h 



h [96] 
[97] 
[70] 
[71] 



+e- - 


h 


+e- - 




+e- - 


h 


+e- - 


h 


+e- - 


h 


+e- - 


h 


+e- - 


h 


+e- - 


h 



[24] 
[24] 
[24] 
[24] 
[58] 
[93] 
[94] 
[95] 



EcM [GeV] 



22 
26 
27.4 
63 



3.1 

5 

4-7 
29 
29-37 
91 
91 
91 



Parameter 



ra [fm] 



0.51 ± 
0.58 ± 
0.54 ± 
0.41 ± 



0.01 
0.07 
0.01 
0.02 



0.53 ± 
0.55 ± 
0.45 ± 
0.64 ± 
0.52 ± 
0.66 ± 
0.58 ± 
0.65 ± 



0.03 
0.03 
0.04 
0.06 
0.07 
0.05 
0.05 
0.07 



r3/r2 



0.61 ± 0.12 
0.58 ± 0.25 
0.45 ± 0.07 
0.50 ± 0.13 



0.65 ± 0.17 
0.65 ± 0.20 
0.63 ± 0.18 
0.77 ± 0.25 
0.59 ± 0.21 
0.80 ± 0.16 
0.73 ± 0.13 
1.00 ± 0.12 



states or they belong to an identical di-baryon state. The most striking feature that emerges 
from the measured data listed in Table 4 is the very small r values obtained for the identical 
di-baryon systems which lie in the vicinity of ~0.15 fm, way below the values obtained for the 
mesons and also much smaller than the proton charge and nuclear radius [105] of ~ 0.83 ±0.02 
fm. Finally the fact that the data in Table 4 utilised the same reaction at the same centre of 
mass energy, affords a unique opportunity to study the emitter dimension r as a function of 
the hadron mass. This feature is dealt with in the following sections. 

6 The r dependence on the hadron mass 

The r2 dependence on the hadron mass, r(m), shown in Fig. 14, exhibits a relatively strong 
decrease in r as the hadron mass increases, a feature which however currently rests mainly 
on the results obtained for the baryon-pairs. The possibility that the origin of this behaviour 
arises from purely kinematic considerations has been explored in Ref. [106] with the conclusion 
that this by itself cannot account for the rather sharp decrease of r2 with m. Thus this r2 
dependence on the hadron mass has to be faced by all non-perturbativc QCD models which 
attempt to describe the hadron production process. In particular r(m) poses a challenge to the 
string approach which also constitutes a basis for the Lund model of hadronisation described in 
details in [47], which expects that dr/dm > in its rudimental form. Whereas the somewhat 
smaller value of r{mK) as compared to that of rim.^), may still be acceptable, the very small 
r2 values extracted from the A and proton pairs cannot apparently be accommodated within 
the Lund string approach [107]. 
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Figure 12: The ALEPH results from the spin-spin FDC analysis of the AA and AA pairs observed in 
the hadronic decays [100] given at several Q regions and presented in terms of dN/dy where y is 
the cosine angle between the two decay protons. The lines represent the best fit of Eq. (28) to the 
data. 

6.1 The r(m) description in terms of the Heisenberg relations 

The maximum of the BEC and FDC effects, from which the dimension r is deduced, occurs 
when the Q value of the two identical hadrons of mass m approaches zero i.e., the hadrons are 
almost at rest in their CMS which also means that the three vector momentum difference Ap 
approaches zero. This motivated the attempt to link the r(m) behaviour to the Heisenberg 
uncertainty principle [102] . From these uncertainty relations one has that 



ApAr = 2iivr 



m vr 



he , 



(40) 



where m and v are the hadron mass and its velocity. Here n is the reduced mass of the di- 
hadron system and r is the geometrical distance between them. The momentum difference Ap 
is measured in GeV and r = Ar is given in Fermi units while he = 0.197 GeV fm. From Eq. 
(40) follows that 

he he 

r = — = — . (41) 
mv p 

Simultaneously one also utilises the uncertainty relation expressed in terms of time and energy 



AEAt ^ —At ^ h , 
m 



(42) 



where the AE is given in GeV and At in seconds. Thus one has 



p 



hm/At so that p 



h m/At . 



(43) 
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Figure 13: Left: The correlation function C{Q) obtained by the ALEPH collaboration [100] using 
Method II for the AA pairs for three different reference samples. A) Monte Carlo, B) and C) two, 
somewhat differently treated, mixed events samples. Right: The OPAL C{Q) distribution of the pp 
system taken from Ref. [101]. The continuous lines represent the best fit of Eq. (29) to the C{Q) 
results. 



Inserting this last expression for p in Eq. (41) one finally obtains 



r{m) 



hc/Jh/At 



zVhAt 



m 



(44) 



If one further assumes that At is independent of the hadron mass and its identity and just 
represents the time scale of strong interactions, of the order of At = 10~^^ seconds, then 
r(m) = A/y/m with A = 0.243 fm GeV^^^. This r dependence on m is shown in Fig. 14 
by the thin solid line. The sensitivity of r(m) to the value of At is represented by the dashed 
lines in the same figure. The upper and the lower dashed lines correspond respectively to the 
values of At = 1.5 x 10"^'* and 0.5 x 10~^^ seconds. A fit of Eq. (44) to the data points plotted 
in Fig. 14 yields At = (1.2 ± 0.3) x lO"^'* seconds where 0.3 represents the statistical error. 
As seen in the figure, the dependence of r on the mass m, as determined from the Heisenberg 
relations, follows within errors the measured values and thus offers a natural explanation for 
the experimental found hierarchy r^^r > ^kk > fpp ~ ''^aa- 



6.2 QCD description of r(m) via the virial theorem 

In the previous section it has been shown that the observed hierarchy dr/dm < has a 
natural and a rather general explanation. Here we would like to demonstrate that such a 
dependence could shed light on the character of the "soft" interaction or, in other words, on 
the non-perturbative QCD which is responsible for the interaction at long distances. In that 
we follow [102] and use the semi-classical approximation as in Sec. 6.1, which means that the 
angular momentum is given by 

£ = \pi - P2\bt = 2pbt ^ he, (45) 
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Table 4: Experimental results for r and A extracted from Bose-Einstein and Fermi-Dirac correlations 
of boson and baryon-pairs present in the hadronic decays produced in e+e~ annihilations in the 
LEPl collider. 



h h 


A2 


r2 [fm] 


Experiment 


^± ^± 

TT TT 




0.78 ±0.01 ±0.16 




\\ 021 




0.82 ± 0.11 ± 0.25 


0.48 ± 0.04 ± 0.07 




r7oi 
[79J 




0.82 ± 0.22 ^nV) 


0.56 ± 0.08 loo6 


OPAL 


f80l 




1.14 ±0.23 ±0.32 


0.76 ±0.10 ±0.11 


OPAL 


[103] 




0.96 ±0.21 ±0.40 


0.65 ±0.07 ±0.15 


ALEPH 


[104] 




0.61 ±0.16 ±0.16 


0.55 ±0.08 ±0.12 


DELPHI 


[79] 


p p 




0.15 ±0.02 ±0.04 


OPAL 


[101] 


A A 




0.11 ±0.02 ±0.01 


ALEPH 


[100] 


A A 


Spin Analysis 


0.19 lg:^?±0.02 


OPAL 


[98] 




a 


0.11 tni±o.oi 


DELPHI 


[99] 




a 


0.17 ±0.13 ±0.04 


ALEPH 


[100] 



where ht is the impact parameter. To estimate the value of p one uses the virial theorem [108] 
which provides a general connection between the average values of the kinetic and the potential 
energies (T and V), namely 

2(Ti) = (6* • \ItV{r) ) , (46) 

where t denotes the transverse direction. Since the motion in the transverse direction is always 
finite, one can safely apply the virial theorem. Substituting in Eq. (46) the kinetic energy by 
its relation to the momentum, Tt = pi /in, we obtain 

{pD = m{n ■ VtV{r) ) , (47) 

where r stands for the distance between the two particles with being equal to 2bt. 

From Eq. (47) follows that the relation between (p^ ) and r depends crucially on the r 
dependence of the potential energy V{r) . Here we take V{r) to be independent of the interacting 
particles (quarks) mass. This assumption is based on the Local Parton Hadron Duality [109,110] 
concept, which states that one can consider the production of hadron as interaction of partons, 
i.e. gluons and quarks, in spite of the unknown mechanism of hadronisation. This approach 
is supported experimentally by the single and double inclusive hadron production in e+ e~ 
annihilation [109,110]. Since the interaction between partons is independent of the produced 
hadron mass it follows that dV{r)/dm = 0. Based on this assumption and using Eq. (45) and 
Eq. (46), one obtains an equation for the typical size of the emission volume, namely 

(f-V V(r) ) « . (48) 

m 

Further one apphes to Eq. (48) the general QCD potential 

V{r) ^ KT - (49) 
3 r 
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Figure 14: The dependence of r2 on the hadron mass obtained from BEC and FDC analyses of 
identical hadron pairs emerging from the decay. For clarity the data points are plotted at slightly 
displaced mass values. The error bars correspond to the statistical and systematic errors added in 
quadrature. The thin solid line represents the expectation from the Heisenberg uncertainty relations 
setting At = W~^^ seconds [102]. The upper and the lower dashed lines correspond to At = 1.5 x 10 
and 0.5 x 10~^^ seconds respectively. The thick solid line represents the dependence of r2 on the hadron 
mass as expected from a QCD potential given by Eq. (49). 

widely used to derive the wave functions and decay constants of hadrons [111, 112]. Setting 
the potential parameters values to k = 0.14 GeV^ = 0.70 GeV/fm , as = 27r/9 ln(5 + 7/r) 
with 6 = 2 and 7 = 1.87 GeV^^ = 0.37 fm as derived from hadron wave functions and decay 
constants [111], the resulting dependence of r on the hadron mass, given by the solid thick 
line in Fig. 14, is seen to follow rather well the the data. Moreover, this QCD expectation 
essentially coincides with the prediction derived from the Heisenberg uncertainty relations (thin 
solid line in the same figure). 

7 Results from multi-dimensional Bose-Einstein correla- 
tion analyses 

The study of multi-dimensional BEC has gained considerable interest in recent years, in par- 
ticular in heavy-ion collisions sector, but also in particle reactions. In the framework of the 
Lund hadronisation model it is expected that the emitter shape will deviate significantly from 
a sphere and that an elongation' of the source should occur. Three LEP experiments, in their 
study of the process e"'"e~ — >• j*/Z^ — > hadrons, have analysed their data in the Longitudinal 
Centre-of-Mass System with the aim of extracting values for the longitudinal (= r^), trans- 
verse rrp and side r<.jf;e dimensions (see Sec. 2.4). These are summarised in Table 5 together 
with the results from the ZEUS experiment at HERA ep collisions which confirm that the ratios 
^side/^z and tt/tz are significantly smaller than one. Missing however is a control experiment 
that will show that in the case of a genuine spherical emitter the multi-dimensional analysis 
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procedure will in fact yield the results that Vgidel'^z ~ fr/fz'^ 1- 

A somewhat surprising result which came out of the BEC multi-dimensional analysis of 
pion-pairs is the dependence of r^, as well as Tside and r^, on the transverse mass rriT- An 
example of this dependence on my is the one reported by the DELPHI collaboration [113] 
which is shown in Fig. 15. This rz{mT) behaviour is very similar to that seen in heavy ion 
collisions where it is often attributed to nuclei collective effects which should be absent in 
particle reactions and in particular in e+e^ annihilations. A theoretical understanding in terms 
of the Lund model of this effect in particle reactions is apparently still lacking even if only for 
the reason that the inclusion of the BEC effect into the model is a highly non-trivial task and 
still seems to need further improvements [114, 115]. 

Table 5: Evidence for the elongation of the source from experimental results for the ratios rside/'''z 
and Tj^/r^ obtained from multi-dimensional BEC analyses of identical charged pions. The first and 
the second errors are respectively the statistical and systematic ones. 



Reaction 



[GeV] 



Experiment 



^ side I "f z 



e+e- ^ Z 
e+e- ^ Z 
e^e" — > Z — >h 



91 
91 
91 



L3 [37] 
OPAL [39] 
DELPHI [38] 



0.80 ± 0.02lHi 
0.82 ± 0.02lg:g^ 



0.62 ± 0.02 ± 0.05 



e p — > e /i 



110 < g2 [GeV^] 



ZEUS 



[16] 



0.67 ±0.08 



7.1 rz;(mT) description in terms of the Heisenberg relations 

We have seen that there is a a great similarity between the dependence of the 1-dimensional r2 
on the hadron mass m and the dependence of the on the transverse mass mx- This is not 
surprising once one realises that the longitudinal distance Ar^, and Ap^, the difference in the 
longitudinal momentum of the two hadrons in the LCMS, are also conjugate observables which 
obey the uncertainty relation [116] 

Ar.Ap, = he . (50) 

Here Ap^ is measured in GeV, Ar^ = is given in Fermi units and he — 0.197 GeV fm. In 
the LCMS one has 

Ap^ = 2fiVz = Pz , 

where = m/2 is the reduced mass of the two identical hadrons of mass m and the longitudinal 
velocity of these hadrons. Thus 

he , , 

Tz = — . 51 
Pz 

Simultaneously one considers the uncertainty relation expressed in terms of energy and time 

AEAt = h , (52) 
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Figure 15: Preliminary results of DELPHI [113] for the dependence of the longitudinal emitter 
dimension r^(7r7r) on the transverse mass mT(7r7r) in hadronic decays. The data are compared with 
the expression for given in Eq. (59) setting At to the best fitted value of 2.1 xlO~^^ sec (continuous 
line) and the expectation for 1.0 xlO~^^ sec (dashed line). 



where AE is given in GeV and At in seconds. In as much as the total energy E of the two- 
hadron system is given by their mass and their kinetic energy, i.e. the potential energy can be 
neglected, one has 

E ^Y.^m'^+Pl+Ply + P'. = EV<t+pL (53) 

i=l i=l 

where rrii T and m2,T are the transverse mass of the first and second hadron. As Qz decreases 
the longitudinal momentum vanishes so that one can, once pi < mj^j,, expand the hadron 
energy E in terms of p^/'^ir retain only the two first terms. 



E = J2^hT, 



Next one can order the identical bosons so that mi,T > Tn2,T and define 

orriT = 77 ris while rriT = — ■ — — ■ (55) 

rrij^ — {dm.T) 2 

Inserting these relations into Eq. (54) one gets, after a few algebraic steps that 

E ^2mr+ ^ ^2 ■ (56) 

rrij' — [omTj 

Thus as long as m\ ^ {dmrY the total energy of the di-boson system is approximately equal 
to 

E K 2mT + — . (57) 
rriT 
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In the 2-dimensional analysis the estimation is deduced from the BEC enhancement in the 
longitudinal direction, keeping rrix fixed and letting and therefore also Pz, approach zero. 
Thus as long as the uncertainty in 2mT is much smaller than that oi pl/rriT one has 

2 

AEAt = ^At = h . (58) 
rriT 

Combining Eqs. (51) and (58) one finally has 



, , cy nl\t . ^. 

rzirnT) ~ — ; . (59) 

This last equation is identical to the one given in Eq. (44) for the dependence of the emitter 
dimension on the boson mass when r and m are replaced by Tz and rriT- A fit of Eq. (59) to the 
data shown in Fig. 15 yields At = (2.1 ±0.4) x lO^^'' sec so that r^(mT) = 0.354/yt[Ge\^ 
fm. This At value is compatible with the corresponding value of (1.2 ± 0.3) x 10~^^ sec ob- 
tained in [102] from a fit to r(m) in view of the considerable spread of the r values obtained 
from the 1-dimensional BEC analyses carried out with the LEPl data. In heavy-ion collisions 
the longitudinal range was also observed to be inversely proportional to the square root of 
rriT [12] namely, 2 / yJniTlGeV] fm where the difference between the proportionality 

factor of 0.354 and 2.0 can be accounted for by the difference in the extent of the heavy ion 
targets as compared to that existing in the e'^e" annihilation. 

The dependence of the transverse dimension tt on the two-pion transverse mass ttit has 
also been measured in the hadronic Z° decays [117]. Here again the transverse range was found 
to decrease as increases. However unlike which is a geometrical quantity, is a mixture 
of the transverse radius and the emission time difference between the two hadrons so that a 
straightforward application of the uncertainty relations is improper. 

7.2 Application of the Bjorken-Gottfried relation to r(mT) 

Another approach to account for the dependence of and on the mass of the hadron 
has been taken in Rcfs. [118, 119]. In that approach the parameters characterising the particle 
emission region are mass independent and the observed change in radii is solely the consequence 
of the momentum-position correlation as expressed in the Bjorken-Gottfried condition [120- 
123] . In this hypothesis there exists a linear relation between the 4-momentum of the produced 
particle and the time-space position at which it is produced, namely 

g^u = Xx/^ which implies A = uit/t , 

where mr is the transverse mass and r = y/t^ — z"^ is the longitudinal proper time after the 
collision. In the framework of this hypothesis relations between the longitudinal and transverse 
radii and the transverse mass can be derived. These are given in Fig. 16 by the dashed bands 
where they are seen to adequately describe the results obtained from the 1-dimensional BEC 
analyses carried out with the LEPl data. 

7.3 r(m) and the interatomic separation in Bose condensates 

The observation that the dimensions r and derived from BEC analyses are seen respectively 
to be inversely proportional to \/m and y^mr, naturally arouses the curiosity whether similar 
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Figure 16: (=^||) and i=R±) as a function of the transverse mass niT taken from Refs. [118,119]. 
The dashed bands were calculated from a model based on the Bj or ken-Gottfried relation. The data 
points are taken from BEC analyses of hadronic decays. 

behaviour does also exist in other phenomena associated to the Bose-Einstein statistics. To this 
end we turn here to the Bose condensate states where we will show that the average interatomic 
separation is inversely proportional to the square-root of the atomic mass. In this presentation 
we follow closely reference [116] and first summarise briefly, for the benefit of the non-expert 
reader, the main features of the Bose condensates and focus on their interatomic separation. 

When bosonic atoms are cooled down below a critical temperature Tb, the atomic wave- 
packets overlap and the equal identity of the particles becomes significant. At this temperature, 
these atoms undergo a quantum mechanical phase transition and form a Bose condensate, a 
coherent cloud of atoms all occupying the same quantum mechanical state. This phenomenon, 
first surmised by A. Einstein in 1924/5, is a consequence of quantum statistics [124,125]. De- 
tailed theoretical aspects of the Bose-Einstein condensation can be found in Ref. [126] and its 
recent experimental situation is described in Ref. [127]. Concise summaries, aimed in particular 
to the non-expert, both of the experimental aspects and the theoretical background, can be 
found in Refs. [128, 129] . To form Bose condensates, extremely dilute gases have to be cooled 
down below the critical temperature Tb, so that the formation time of molecules and clusters 
in three-body collisions is slowed down to seconds or even minutes to prevent the creation of 
the more familiar transitions into liquid or even solid states. 

The existence of Bose-Einstein condensation was first demonstrated in 1995 by three groups 
[130-132] in cooling down rubidium, sodium and lithium. Typical temperatures where Bose 
condensates occur are in the range of 500 nK to 2 /xK with atom densities between 10^^ and 
10^^ cm^^. The largest sodium condensate has about 20 million atoms whereas hydrogen con- 
densate can reach even one billion atoms. 

Let us now consider a dilute homogeneous ideal gas of N identical bosonic atoms of spin 
zero confined in a volume V. These atoms occupy energy levels, handled here as a continuous 
variable, which are distributed according to the Bose-Einstein statistics. We further set the 
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ground state to be Sq — 0. If A^o is the number of atoms in this ground state and N^.^ is the 
number of atoms in the excited states then = A'^o + N^^. For a homogeneous ideal gas of 
identical bosonic atoms it can be shown (see e.g. Ref. [133]) that at a low temperature T one 
has 

= 2.612 V (60) 

where V is the volume occupied by the atoms of mass m. Since Tg is defined as the temperature 
where almost all bosons are still in excited states, we can, to a good approximation, equate N 
with Nex- That is 



AT = 2.612 . (61) 

For T < Tb one obtains from Eqs. (60) and (61) that the number of atoms A"o which are in 
the condensate state is, 

rp . 3/2" 



1 - 



Tr 



■J 



(62) 



where T is the temperature of the atoms lying at the excited energy states above the condensate 
energy level eo = 0. The atomic density of the Bose gas at very low temperatures, T/Tb « 1 
where N ^ No, is then given by 

where k is the Boltzmann constant and p, the atomic density, has the dimension of L~^. From 
this follows that p~^^^ is the average interatomic separation in the Bose condensate. At the 
same time the thermal de Broglie wave length is equal to 

Combining Eqs. (63) and (64) one has for the state of a Bose condensate the relation 

pXIb ~ 2.612 . (65) 

Thus the average interatomic distance in a Bose condensate, (Ibe-i is equal to 

dBE = P"'/' ~ AdB/1.378 . (66) 

Next one can consider two different bosonic gases, having atoms with masses mi and m2, 
which are cooled down to the same very low temperature Tq, below the critical temperature Tb 
of both of them. In this case one produces two Bose condensates with interatomic distances 



From this follows that when two condensates are at the same fixed temperature Tq one has 



dBE{m2) 




(68) 
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which is also the expectation of Eq. (44) for the r dependence on the mass of the hadron pro- 
duced in high energy reactions provided At is fixed. Finally it is interesting to note that in as 
much as one is justified, since almost all the atoms occupy the same quantum mechanical state, 
to replace in Eq. (67), at very low temperatures, kTo by AE and use the uncertainty relation 
AE — h/At one derives the expression for r(m) as given by Eq. (44) multiplied by the factor 
-\/27r/1.378. Thus in Bose-Einstein condensates r(m) measures the interatomic separation and 
not the dimension of the condensate ensemble, a deduction which will be further referred to in 
Sec. 8. 

In relating the condensates to the production of hadrons in high energy reactions one should 
keep in mind that the interatomic separation proportionality to 1/ ■sjmatom does not necessarily 
imply that it should also apply to hadrons produced in high energy reactions. Common to 
both systems is their bosonic nature which allows all hadrons (atoms) to occupy the same 
lowest energy state. In addition the condensates are taken to be in a thermal equilibrium state. 
Among the various models proposed for the production of hadrons in high energy e+e" and 
nucleon-nucleon reactions some attempts have also been made to explore the application of a 
statistical thermal- like model [134]. Whether this approach will eventually prevail is at present 
questionable. Finally the condensates are taken to be in a coherent state. In the case of hadrons 
one is able to measure r via BEC only if the chaoticity factor A in Eq. (9) is different from 
zero i.e., only if the source is not 100% coherent. However so far there is no evidence from 
BEC studies in hadron production for a dependence of r on A apart from that introduced by 
the correlated errors between Ar and AA produced by fitting Eq. (9) to the data. 

8 On the relation between r and the emitter size 

As it has been pointed out in the introduction, the HBT interferometry in astronomy and its 
equivalent, the GGLP correlations in its application to heavy ion collisions, have been rather 
successful in estimating respectively the dimensions and configurations of stellar objects and 
the geometrical extent of nuclei. In the elementary particle sector it looks as if the r dimension 
is very little if at all dependent on the type and energy of the reaction. Furthermore e^e~ 
annihilations, there seem to be several features which pose a challenge to the general belief 
that r measures a quantity which represents the emitter volume and corresponds to its radius 
in the spherical approximation. 

In the Lund string model one expects that the extension of the emitter of hadrons in ele- 
mentary particle reactions will increase with multiplicity and hence with energy and certainly 
in e^e^ annihilation at 91 GeV will spread over more than several fm. From the Bose-Einstein 
correlation analyses there is no evidence for such a behaviour. In fact, the r values obtained 
from T^^i:^ pairs produced in e+e~ annihilation seem to be rather independent of energy having 
values in the vicinity of 0.75 fm (see Table 1). 

Moreover, recent experimental results show that the dimension r, measured in e^e~ anni- 
hilations on the Z° mass energy, does depend on the hadron mass. As shown in Sec. 6, the 
r values at the proton and A masses decrease to the very low value of ~0.15 fm as compared 
to approximately 0.75 fm for pions and kaons. This behaviour of r(m) is very similar to the 
interatomic separation dependence on the atomic mass in equal temperature Bose condensates. 
In as much as one is justified to infer from the Bose condensates features to the hadron pro- 
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duction sector, r(m) apparently does not measure the emitter radius but rather the distance 
between two neighbouring identical hadrons. 



In any case, the experimental results found for the dependence of r on the hadron mass are 
not reconcilable with the notion of a unique emitter radius for all outgoing hadrons. However 
the possibility that baryons are emitted from a much smaller source than pions and kaons has 
its own difficulty [135] if one considers the source energy density e. To illustrate this problem 
let us assume that the emitter source is a sphere of radius r, as determined from the correlation 
analysis, filled homogeneously with at least the energy which corresponds to the sum of the 
rest masses of the two outgoing identical hadrons. Then, for a given hadron mass rrih and its 
emitter radius r^, the energy density of the source is given by 



3 2mh . . 

^measured . q V^''/ 

47r rl 



which is plotted in Fig. 17 for the pion, kaon, proton and A pairs. The dashed lines in the 
figure are the predicted dependence of the energy density on the hadron mass computed from 
Eq. (70), which is based on the expression of r(m) given in Eq. 44, namely 

_ 3 m^' 

Whereas the energy density of the pion and kaon pairs are still lying within an acceptable range, 
the energy density of the proton and A emitter sources reach the very high value of about 100 
GeV/fm3. 
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Figure 17: Energy density of the hadron emitter, calculated from Eq. (69), as a function of the 
hadron mass assuming a spherical emitter having the radius r which was extracted from the BEC 
analysis. The dashed lines are the expectations from Eq. (70) setting At to plus and minus three 
standard deviations away from its median value of 1.2 x 10^^"^ seconds [135]. 

In trying to cope with these observations, it may be helpful to try and identify the un- 
derlying variable or variables, apart from the hadron mass, which r depends on. A partial 
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answer to this quest [136] , which is briefly outhned in the foUowing section, is provided by the 
relation between the BEC analysis and the genuine particle correlation method which utilises 
the factorial cumulant moments (see Sec. 2.3). 



8.1 The r dependence on the number of emitting sources 

Let us consider a simplified description of heavy ions (AA) interactions at high energies which 
is given in terms of several distinct, say S, primary nuclcon-nucleon (nn) collisions leading to 
hadrons. In this approach the hadrons produced in a single primary nucleon-nucleon reaction 
are emitted from an S'„„ source. The hadron produced in heavy ion (AA) collisions are then 
coming from Saa sources so that Saa/ Snn—^- Since here the AA collisions are described in 
terms of the single primary nucleon-nucleon reactions one can define, without loss of generality, 
Snn as a unit source so that Saa = S. 

Next the tt^tt^ pairs emerging from a single primary nucleon-nucleon reaction and those 
emerging from the heavy ion collisions are subject to BEC analyses and the correlation functions 
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Ct\Q) = 1 + As=ie-« '■1=1 and C^{Q) = 1 + A^e-^^ ^: 

are fitted, each to its corresponding data distribution. In general for every 5* vahic the resulting 
rs can be extracted from the fitted correlation functions by integration C2 — 1 over Q from 
zero to infinity so that 
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CUQ)-l\dQ=l Xse-^'^sdQ =1^, (71) 
■'Jo z rs 

where A5 is, as customary in the BEC analyses, taken to be independent of Q. 

For the two-particle correlation there exists a relation (see e.g. Ref. [137]) between the 
cumulant K2 and the BEC function C2, namely that 

C2(Q) = 1 + Xe-"-"^" = 1 + ^2 . (72) 

In addition it has been shown [138-140] that for two particle correlation the cumulant of S 
sources is, to a very good approximation, related to the value of the corresponding cumulant 
of one source, -fi'f"^, by 

Thus the BEC function Cf is equal to 

By integrating Eq. (73) over Q from zero to infinity one finally obtains that 

rs = S^rs=i , (74) 

where S like the A parameters is taken to be independent of Q. Consequently if the increase 
of an emitter volume is coupled to an increase in the number of boson sources, then the BEC 
extracted rs will grow as well. Further to note is that if in heavy- ion collisions SXs/Xs=i is 
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proportional to A^^^ then r will be proportional to the nucleus radius. As for the HBT analysis 
applied in astronomy, the proportionality between the stellar object dimension and the number 
of photon sources is self-evident. 

Prom the discussion above it is obvious that the choice of the unit source determines the 
predicted dependence of r on the various features, such as type and energy, of the hadron (or 
photon) emitting reaction. In the case of heavy ion collisions the choice of the unit source to 
be a single primary nucleon-nuclcon collision is a natural one. In the case of e+e^ annihilation 
one may consider the possibility that S^e itself is the unit source. Under this assumption the 
number of sources does not change with energy and r will remain constant. This prediction 
seems to be consistent with the measured r values allowing for the fact that the various ex- 
periments used different analysis methods and reference samples (see Sec. 5.1). On the other 
hand an increase of r as a function of the hadron multiplicity has been observed in the BEG 
analysis [22] of equally non-zero charged pion-pairs emerging from the decay. Furthermore 
the same experiment also found that the extracted r value from the Z° decays into three hadron 
jets was by about 10% higher than that obtained from the Z° decay events having only two 
hadron jets. Similar rise of r with the average hadron multiplicity has also been reported by 
several hadron-hadron collision experiments carried out at centre of mass energies above 30 
GeV [73,141-143]. 

These last findings indicate that apparently for the e'^e" annihilations a reasonable choice 
for the unit source is the single hadron jet. The shortcoming of such a choice is however the 
fact that the charge multiplicity of a hadron jet depends on its energy and on its origin, i.e. if it 
is a quark or a gluon [144-146] , which currently impedes a detailed quantitative prediction for 
the behaviour of r as a function of the centre of mass energy and the multiplicity. Regardless 
of these hindrances, the adaptation of the hadron jet as a unit source provides nevertheless a 
qualitative explanation for the increase of r with multiplicity and also expects in e'^e~ annihi- 
lations a moderate dependence of r on the e+e" centre of mass energy. 



9 Summary and conclusions 

The 1-dimensional Bose-Einstein correlation enhancement measured in identical pion-pairs close 
in time-phase space, is found to be a universal phenomenon present in a large variety of parti- 
cle reactions and over a wide range of energies. Nevertheless, there are still missing systematic 
BEG studies concentrating on the same reaction at different energies, and/or different reactions 
at the same energy, which will allow a comprehensive investigation of the dependence of r on 
the variables that characterise the particle interactions. Notwithstanding the fact that these 
systematic studies are still missing and that BEG analyses cannot currently be considered as 
precision measurements, the extracted r values, which are taken to represent the hadron emit- 
ter source, seem to be very little, if at all, dependent on the reaction type or on its energy and 
they are seen to fluctuate in the range of about 0.6 to 1.0 fm. 

The essentially constant r value seems to be in contradiction with the results obtained 
from the BEG studies of heavy- ion collisions and from the HBT analyses of stellar objects. 
An insight into the origin of this difference is offered by the relation between the BEG and 
the cumulant moments analysis where the introduction of the so called hadron source plays a 
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central role. Prom this relation follows that the value of r is proportional to the number of 

active hadron sources which participate in the reaction. In particle reactions the proposal to 
equate the number of sources with the number of hadron jets offers an understanding of the 
r dependence on multiplicity. At the same time however the difficulty in defining a standard 
hadron jet prohibits presently any quantitative results. 

The genuine higher order BEC analyses are presently restricted to three identical charged 
pions. Assuming a Gaussian hadron emitter, the expected relation r2/\/3 < ra < ra/v^ be- 
tween the two and three-pion BEC dimension parameters, is found to be consistent within 
errors with the measurements carried out in e'^e~ annihilation and hadron collisions. The 
direct measurement of the genuine correlations of four and more pion systems has not been 
realised so far. On the other hand the extracted r values from the , so to speak, non-genuine 
correlations, are model dependent and as such less reliable. 

Recently first results from 2 and 3-dimensional BEC analyses of pion-pairs have been re- 
ported. These verified the Lund string model expectation that the hadron emitter source 
deviates from a sphere and is elongated in the thrust (or sphericity) direction of the multi- 
hadron final state. In the approximation that the emitter has an elongated prolate ellipsoid 
shape, where the long axis is in the thrust direction and rr is the axis perpendicular to it, 
the ratio r^/r^ is found be in the range of ~0.7 to ~0.8. 

The extension of the BEC analysis to the fermion sector was shown to be feasible via the 
Pauli exclusion principle. Pirst results from the so called Permi-Dirac correlation analysis of 
the AA and proton— proton pairs, emerging from the hadronic Z° decays, have demonstrated 
that the dimension r is strongly dependent on the mass of the identical hadrons. It has been 
further shown that this dependence can be accounted for by the Heisenberg uncertainty rela- 
tions, or alternatively by QCD derived potential, to yield the relation r{m) oc l/^/m. Similar 
relation is seen to be valid also in the 2-dimensional BEC analysis when r is replaced by the 
and m by the transverse mass mj. It is interesting further to note that the interatomic sepa- 
ration in the Bose condensates having the same low temperature, is proportional to 1/ ^Jraatom ■ 

So far the PDC analyses were restricted to the hadronic decays, therefore it is essential 
that they will be repeated in other particle reactions to verify that indeed r of the A and the 
proton pairs is in the vicinity of 0.15 fm. This small r value, derived from the PDC analyses 
of baryon-pairs, poses a challenge to the concept that r represents the hadron emitter radius. 
If one does insist on keeping the radius emitter interpretation for r then one encounters the 
situation of an unreasonable high energy density value, of some 100 GeV/fm^, of the baryons 
source. Thus the very small taa and r^^ values indicate that our level of understanding the 
production of hadrons, and in particular baryons, is apparently far from being satisfactory. 
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